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Abstract 

This paper concerns the reconstruction of possibly complex-valued coefficients 
in a second-order scalar elliptic equation posed on a bounded domain from knowl- 
edge of several solutions of that equation. We show that for a sufficiently large 
number of solutions and for an open set of corresponding boundary conditions, all 
coefficients can be uniquely and stably reconstructed up to a well characterized 
gauge transformation. We also show that in some specific situations, a minimum 
number of such available solutions equal to /„ = ^n{n + 3) is sufficient to uniquely 
and globally reconstruct the unknown coefficients. This theory finds applications 
in several coupled-physics medical imaging modalities including photo-acoustic 
tomography, transient elastography, and magnetic resonance elastography. 

1 Introduction 

We consider the general second-order elliptic equation: 

V ■ aVu + b-Vu + cu = in X, u = f on dX, (1) 

with complex-valued coefficients such that a is a symmetric, possibly complex-valued, 
tensor verifying the ellipticity condition for > 0: 

ao|er<e-5Rae<ao'ler (2) 

and with c such that the above equation admits a unique solution in H^{X) for / G 
H2{dX). Here X is an open bounded domain in M" with smooth boundary dX. We 
assume throughout this paper that a, b, c and V ■ a are of class C^'^^^X) for some a > 0. 
Elliptic regularity results [TU Theorem 6.3.7] then ensure that the solution u is a. strong 
solution of class C^'"(X) when the boundary condition is of class C'^'"'{dX). 



*Dcpartment of Applied Physics and Applied Mathematics, Columbia University, New York NY, 
10027; gb2030@columbia.edu 

'''Department of Mathematics, University of Washington, Seattle, WA, 98195 and University of 
California, bvine, CA, 92697; gunther@math.washington.edu 



1 



We assume that we have access to internal functionals given by the complex valued 
solutions Uj{x) of ([1]) for a given set of boundary conditions / = fj for 1 < i < /• 
The main objective of this paper is to show that when / is sufficiently large, then the 
coefficients (a, b, c) can be uniquely and stably reconstructed from knowledge of the 
corresponding (iij)i<j</ up to a natural gauge transformation. 

In some specific situations, for instance when (a, 6, c) is close to (ao,0,0) for uq a 
constant complex- valued matrix satisfying ([2]), then the reconstructions can be achieved 
for I = In = |n(n + 3) coefficients, which corresponds to the (complex) dimension of 
the unknown coefficients (a, b, c) up to the (one parameter) gauge transformation. 

The mathematical description of the measurement operator considered here and the 
main results of the paper are presented in section O The proofs of the main results are 
detailed in sections [3] and IH 

The theory presented in this paper finds some applications in several recent coupled- 
physics (also called hybrid or multi-wave) imaging modalities that are extensively stud- 
ied in the bio-engineering community. These methods arose as an attempt to devise new 
imaging modalities that combine high resolution with high contrast. In section [5l we 
consider three such modalities: Photo-Acoustic Tomography (PAT), Transient Elastog- 
raphy (TE) and Magnetic Resonance Elastography (MRE). In PAT and TE, the high 
resolution mechanism is ultrasound. In MRE, it is Magnetic Resonance Imaging. Our 
results show that in these three imaging modalities, all the medium parameters (some 
of which displaying the high contrast we are after) can be reconstructed uniquely and 
stably. The reason is that these coefficients have a structure that allows us to uniquely 
reconstruct the gauge as well. In PAT, the main novelty of our result is that we can re- 
construct general anisotropic diffusion coefficients. In a scalar model for TE and MRE, 
the main novelty is that we can not only reconstruct anisotropic coefficients but also 
complex- valued coefficients that account for possible dispersion and attenuation effects. 

For the application to quantitative PAT, we refer the reader to e.g. [SI El [7] and the 
references there. For applications in TE and MRE, we refer the reader to e.g. to [12] 
and its list of references. For general references to the theory of hybrid inverse problems 
devoted to the mathematical analysis of similar coupled-physics imaging modalities, we 
refer the reader to e.g. [21 HI [TB] . 

General real-valued anisotropic tensors were reconstructed in the context of ultra- 
sound modulation in dimension n = 2 in [13]. Note that the reconstruction of anisotropic 
coefficients a from boundary value measurements of u (as in, e.g., the Calderon problem) 
can be performed only up to a very large class of changes of variables [20]. Moreover, 
the corresponding stability estimates are of logarithmic type, which corresponds to po- 
tentially drastic amplifications of measurement noise during the reconstruction. The 
results of this paper show that the availability of internal functionals such as those in 
PAT, TE, and MRE, allows one to reconstruct a larger class of coefficients and with 
significantly better stability estimates. 

2 Main results 

Gauge transform. The elliptic equation ([1]) may be recast as 

a : + (V ■ a + 6) • Vtt + c« = 0. (3) 
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Upon multiplying through by the scalar function r 7^ 0, we get 



ttr : V®^u + (V ■ + K) ■ Vu + CrU = 0, 



(4) 



with 



ttr = Ta, br = rh — aVr, Cr = 



rc. 



(5) 



This shows that the coefficients (a, fo, c) can be reconstructed at most up to the above 
gauge transformation. We say that (a, fo, c) ~ (a,-, 6,-, c,-) belong to the same class of 
equivalence if there exists a non vanishing (sufficiently smooth with smoothness de- 
pending on context) function r such that ^ holds. Then we say that (a, 6, c) e c the 
class of equivalence. 

Measurement operator. For / G H^(dX), we obtain a solution u E H^iX) and we 
can define the solution operator 



Note that the solution operator is independent of the element (a, b, c) G c. 

The main measurement operator we consider in this paper is defined as follows. Let 



Here, X is a subset of a Banach space in which the unknown diffusion tensor is defined. 
That space will depend on the context. Also 2) is a subset of H^{X) where the solutions 
to ([I]) are defined. The main objective of this paper is to consider settings in which 9?lf (c) 
for an appropriate choice of f uniquely and stably determines c. 

Main results. The main result of this paper are summarized in the following Theo- 
rem: 

Theorem 2.1 Let c and c be two classes of coefficients with elements {a,b,c) of class 
C"^'"{X) for a > and V ■ a also of class C™'"(X) for m = or m = 1. We assume 
that ([1]) is well posed for the coefficients c. 

Then for I sufficiently large and for an open set (for instance in the topology of 
C^'"(X) ) of boundary conditions f = {fj)i<j<i, then 9?tf(c) uniquely determines c. More- 
over, for any other c as given above, we have the stability results 



for m = 0,1. 

We say that c is in the vicinity of c if {a,b,c,'V ■ a) is in the C°'°'(X) vicinity of 
(a, 6, c, V ■ a) for some elements (a, b,c) E c and (a, b, c) G c. 



H^dX) H\X) 

f ^ u = ej. 



(6) 




||(a, 6 + V ■ a, c) - (a, 6 + V ■ a, c)||i^™,^(x) < C\\Tt^{c) - Tl^{c)\\w^+^.°°{x), 

\\b-b\\Lo^(x) < C\\mf{c)-m^{i)\\w-^,oo^x), 
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Let us assume that c in the vicinity of either: 
(i) (ao,0,0) for some constant diffusion tensor Uq; 

(a) (7(x)/(i„, 0, c(x)) for some scalar coefficients 7 G if2-+^+'^(X) and c E H^^'^^''[X) 
for e > 0; 

(Hi) (7(x), 0, c(x)) for an arbitrary diffusion tensor of class H^^'^{X) and c of class 
H^^^{X) in dimension n = 2. 

Then for / = = ^n{n + 3) and an open set of boundary conditions f, we have that 
njlf(c) uniquely determines c. Moreover, ([8]) holds. 

The proof of the different results described in the above theorem is presented in detail 
in sections [3] and m 

Reconstruction of the gauge. In some situations the gauge in c can be uniquely 
and stably determined. Let us for instance consider the specific, practically important, 
case of coefficients (a, b, c) = (7, 0, a). Then we have the following result: 

Corollary 2.2 Under the hypotheses of the preceding theorem, and in the setting where 
b = 0, we have that 3Jtf(c) uniquely determines (7,0,0). Let us define 7 = rM° where 
M° has a determinant equal to 1 . Then we have the following stability result 

||r- f|Ui.^(x) + ||(MO,c) - (M°,g)|U-(x) < C\\Tlfic)-m^Cc)\\w2,^^x)- (9) 

The proof of the corollary from the results stated in Theorem 12.11 may be found in 
section 13.21 

3 Local reconstruction 

3.1 Reconstruction of equivalence classes 

Let n be the spatial dimension and define In = \n{n + 3). We wish to reconstruct 
an element in the class c = {ar,br,Cr) for r an arbitrary non- vanishing function from 
knowledge of for 1 < i < /„. We perform the reconstruction locally in the vicinity 
of a point Xq. We assume that we have constructed a solution ui such that w.l.o.g. 
Ui{xq) = 1 and by continuity -Ui 7^ in the vicinity of Xq. We then define 

Vj = ^^^^ , 1 < J < -^n — I5 01 = u\a, j3 = V ■ {au\) + u\b 

and find after some algebra that 

a : V'^'^Vj + /3 ■ Vvj = 0, (10) 

in the vicinity of Xq. 

We assume that (Vwi, . . . , Vt'n) forms a frame in the vicinity of Xq. Then H = 
{Hij)i<i,j<n with Hij = Vvi ■ Vvj is a boundedly invertible symmetric matrix with H^^ 
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the coefficients of if ^. Let us define M„ := ^n{n + 1) — 1. Note that = 1 + n + Mn- 
We then write for 1 < m < 

- VVm+n = Q^VVj, = -W'VVm+n " V^;^. (11) 

Here, we use the convention of summation over repeated indices. We then define 
(Of l<j<n 

^f=< 1 j = m + n , l<j</n-l, l<m<Mn. (12) 
l_ otherwise 

Note that we have constructed the complex-valued coefficients in such a way that 

In-l 

J2 OfVvj = 0, l<m<Mn. (13) 
We next construct the symmetric matrices 

= ^ efV^^j, 1 < m < M„. (14) 

j=i 

We assume that the matrices are linearly independent. Since the (complex) di- 
mension of symmetric matrices equals M„ + 1, this implies the existence of a unique 
symmetric, complex-valued, matrix M° such that 

M° : M"^ = Tr(M°M"*) = Som, < m < M„. (15) 

In other words, there exists a unique normalized matrix (M°)* that is orthogonal to the 
constructed M™ matrices for the inner product {A, B) = Ti{A*B). The construction of 
can be obtained as follows. In the constructions presented later in the paper, the 
identity matrix Idn is not in the span of the matrices M"*. We can then use the Gram 
Schmidt procedure to orthonogonahze {Idn, (^'")i<m<A/„) and this way construct the 
matrix M". 

Multiplying f[T(}]) by 6*™" and summing over j yields the M„ constraints 

a : M"' = 0, 1 < m < M„. 

This shows that a is proportional to the now known matrix M°(a;). Note that since the 
real part of a is positive definite, the matrix a cannot be orthogonal to Idn for the inner 
product (■, ■). This justifies the fact that {Idn, (^™)i<m<M„) form a free family. Let us 
define 

^(x) := M°(x). (16) 

Since the matrix a can be reconstructed up to the gauge transformation corresponding 
to the multiplication by an arbitrary complex-valued function, we choose the gauge as 
prescribed above. 

The vector field /3 is then uniquely determined using (fTOll by the explicit formula: 

P = -W^{M'^ : V®^Vj)Vvi. (17) 
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Once [a, (3) are known (up to the gauge transformation), we define 

a=^a, b = ^{(3 — V ■ a), c = (V ■ aViii + 6 ■ Vui). 



(18) 



This provides an exphcit reconstruction of (a, b, c) G c, the class of equivalence, which 
is therefore now known. The above results may be recast as follows 



a=—a, V-a + 6= — (/3-aVMi), c= ( (V ■ a + 6) ■ Viii + a : V^^tii ) . (19) 



We summarize the above results in the following lemma. 

Lemma 3.1 Let Ui for 1 < i < In be solutions of the elliptic equation with boundary 
conditions Ui = fi on dX . Let us define Vi = u^'^UiJ^i for 1 < i < In — ^ and assume 
that: 

(i) in the vicinity Xq of a point xq, we have that Ui ^ on X and Ui{xq) = 1, 

(a) the vectors (Vvi, . . . , Vvn) form a frame on Xq with a matrix Hij = Vvi ■ Vvj with 
uniform bounded inverse on Xq; 

(Hi) the matrices M"^ for 1 < m < Mn constructed above in (fTT |) -( fT^ are linearly 
independent. 

Then the class of gauge equivalence c is uniquely determined by Ttf{c) = (wi)i<i</„ on 
Xq in the sense that {ui)i = {ui)i implies that c = c where Ui are the internal functionals 
obtained by replacing one element in c by one element in c. 

Moreover the reconstruction is stable in the sense that 



\\{a,b + V ■a,c)-{a,b + V ■a,c)\\L^^Xo) < C\\m^{c) - m^{~c)\\w^,oo^Xo), (20) 



for some (a, b,c) E c and (a, b, c) G c. 

The stability result is clear by inspection of the reconstruction procedure: two derivatives 
on V are taken in the reconstruction of the matrices M™' and hence of M° for instance 
by the Gram-Schmidt procedure, a multi-linear operation that preserves errors in the 
uniform norm. The same loss of derivatives is observed in the reconstruction of (a,/3), 
and hence in a, V ■ a + 6 and c as can be seen in ( fT9l) . Note that we similarly reconstruct 
the above coefficients in the W'^'°° sense, and hence b in the L°° sense, when errors are 
bounded in the W^'°° sense as indicated in Theorem 12.11 

Remark 3.2 Note that the procedure described in Lemma lSTll allows one to reconstruct 
In = \n{n + 3) complex-valued coefficients from In complex-valued internal functionals 
of the form Uj above or to reconstruct In real-valued coefficients from In real-valued 
functionals of the same form. 

3.2 Reconstruction of the gauge 

Let us assume that we can reconstruct {oq, b^, cq) G c on a domain X, which may be a 
part of the domain such as Xq above or the whole domain if local reconstructions of c 
are patched together to provide a global reconstruction. Let us then define 
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{ttr, br, Or) = {tOq, T^q - AoVt, TCq) G C 
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an other element of the equivalence class. In this section, we show that prior information 
about the unknown coefficient allows us to uniquely determine the unknown gauge r. 
Note that another way to represent the gauge transform is to realize that 

{ar,a~^br,Cr) = (rao,ao^6o - Vlnr, rco). 

In other words, a~^br — a^^bo is independent of the element in the class c. If a~^br is 
seen as the 1-form {a~^)ijbjdxi, then d{a~^br) is independent of the element in c and 
characterizes the class of equivalence on a simply connected domain X. 
Let us assume that V ■ {a~^br) = $ is known. Then we observe that 

-Alnr = $ - V ■ (tto^feo), 

so that r is uniquely defined on a simply connected domain X if it is known at the 
boundary dX. 

As another practical assumption to reconstruct the gauge, let us assume that b = br 
is divergence free so that V ■ (r^o — «oVr) = 0, or in other words 

-V ■ (aoVr) + V ■ (rbo) = 0. 

Note that uq = u^'^M^. This is an elliptic equation. Provided that all coefficients are 
real-valued and that Oq is uniformly elliptic, then this equation admits a unique solution 
for r when r is known on dX. This is a consequence of the maximum principle [1] that 
does not apply in the case of complex-valued coefficients. 

Let us assume the stronger constraint that b = br = 0- This corresponds to $ = 
above. In fact, this provides the redundant system of transport equations for r: 

Vr = ^6oT, or equivalently V\nT = aQ^bo, 

which admits a unique solution provided that r is known at one point (and admits a 
solution provided that d^aQ^bo) = for a^^bo seen as a 1-form). This reconstruction 
applies for arbitrary complex valued coefficients (a,6, c). In the PAT, TE, and MRE 
applications considered in section [5l the natural setting is with 6 = so that the gauge 
can indeed be reconstructed. 

Consider the specific example of 

V ■ 7Vn = in X, u = f on dX. 

Let us define 7 = rM°. The reconstruction of M° requires taking two derivatives of the 
data. The above equation for r is in fact a redundant system of first-order equations for 
In r in which we gain back one derivative. We thus obtain the unique reconstruction of 
7 with the stability estimate 

\\r - f |Ui,oo(^) + llAfO - M°|Uoo(x) < C||mtf(c) - Mf(c)||w^2,oo(^). 

This proves Corollary 12.21 The reconstruction of the anisotropy is less stable than that 
of the isotropic component. This is consistent with similar results obtained for the 
ultrasound modulation problem; see [13]. 
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3.3 Variations of the coefficients 

The above reconstruction shows that reconstructions are stable with respect to fluctu- 
ations in the measurements 9Jt(c) when the hypotheses of Lemma [3TT] are satisfied. We 
show that such hypotheses are stable with respect to small changes in the parameters 
c. 

Lemma 3.3 Let u and u he solutions of 

V ■ aVu + h ■ Vu + CM = V ■ aVu + h ■ Vu + cu = Q on Xq, 
with Dirichlet conditions u = u = f for f of class C""'"(5Xo). Then 

\\u - n||c"»+2.-(Xo) < C'IKc, V ■ a) - (c, V ■ a)||c™."(Xo), (21) 
for some positive constant C independent of c and oft for c hounded hy M in C'"'". 
Proof. Let w = u — u. We find 

V ■ (a - a)Vu + {b - b) ■ Vu + {c - c)u = V ■ aVw + b ■ Vw + cw. 

The proof then follows from standard regularity results; see [9] for the case of real-valued 
coefficients and [TH Theorem 6.3.7] for the case of complex-valued coefficients. □ 
With m = 0, we deduce that when the hypotheses of Lemma 13.11 are verified for the 
coefficients c, then they are verified with the same boundary conditions for all coefficients 
c that are sufficiently close to c in the sense given above. 

4 Global reconstruction 

We now consider several settings in which global reconstructions of c are possible. We 
refer to the preceding section and the section on applications for reconstructions of the 
gauge under additional information. 

4.1 Global reconstructions close to constant tensor 

We start with reconstructions in the vicinity of a = 6 = 0, and c = 0, where Idn is 
the identity matrix in dimension n. The main interest of this result is that the boundary 
conditions fi are explicit and particularly simple. Moreover, in the case where the 
coefficients (a, fo, c) are complex-valued, the following result shows that the /„ = \n{n+?)) 
unknown complex-valued coefficients may uniquely and stably be reconstructed from 
exactly /„ complex-valued internal functionals. The same proof shows that when all 
coefficients are real-valued, then the /„ = |n(n-|-3) unknown real-valued coefficients may 
uniquely and stably be reconstructed from exactly real-valued internal functionals; 
see remark [3^21 

Theorem 4.1 Let X he the space of c sufficiently close to {ldn.,0,0) in the sense of 
Lemma WTR with m = 0. Let In = ^n{n+3) and fi he the traces on dX of In homogeneous 
harmonic polynomials of degree equal to 1 or 2 (see the proof for the description of the 
polynomials). Then QJtf from X to 2)'^" is injective. Moreover, the stability result fl20|) 
holds for Xq = X . 
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This is therefore a global (in space) stability result but only for classes c that admit an 
element sufficiently close to (/d„,0,0). In fact, the reconstruction works for c close to 
(ao,0,0) for oq arbitrary elliptic as we observe in Theorem 14.21 below. 

Proof. Let iti = 1 be the constant solution. Let then vj = Xj for 1 < j < n 
be the homogeneous polynomials of degree 1. Finally, let us denote by Vij = XiXj 
and Wi{x) = \{xf — xfj^^) the homogeneous harmonic polynomials of degree two for 
^ ^ i < j ^ n. The other such polynomials, such as for instance — Sg) can be 
constructed by linear combination of the polynomials Wi. We have thus constructed 
1 + n + ^n{n — 1) + n — 1 = ^n{n + 3) = /„ harmonic homogeneous polynomials of 
degree less than or equal to 2. 

We verify that 

Vmi = 0, Vvi = ei, Vvij = Xjei + XiCj, ViOi = XjCi - Xi+ie,+i. 

Moreover, 

V^'^Vij = Cj (g) Cj + Cj (g) Cj, V®'^Wi = Ci® ei~ Cj+i ® Cj+i. 
Let us define for 1 < i < j < n and 1 <k <n such that 

{—Xj when k = i 
-Xi when k = j 
otherwise. 

Then we find that the corresponding matrices are defined by 

j^ij ^ ^®'^Vij = ei® tj + ® et. 
Let us now define 0^ for 1 < i < n — 1 and 1 < k < n such that 

{—Xi when k = i 
Xj+i when k = i + 1 
otherwise. 

Then we find that the corresponding matrices M™ are defined by 

M' = V^^Wi = Ci ® - ei+i ® Ci+i. 

The matrices M*-' for 1 < i < j < n and M* for 1 < i < n — 1 form a free family of 
dimension — 1) + n — 1 = |ri(n + 1) — 1 = M„ as can easily be verified. They are 
orthogonal to the matrix M° = Idn- 

Let now c be close to {Idn, 0, 0) in the sense given in Lemma [373] and let the functions 
Ui for 1 < i < /n be solutions of the equation ([1]) with boundary conditions fi that are 
the traces of the harmonic polynomials constructed above. (This means that /i = 1, 
/2 = xi, and so on on dX.) 

By continuity of the solution to ([1]) stated in Lemma 13. 3[ the linear independence 
of the vectors Vvj for 1 < j < n still holds. The linear combinations 6J^ in f|T2|) and 
the matrices M*" in f|T^ constructed by continuity from the case c = [Idn, 0, 0) still 
satisfy (|T3l) and the fact that the matrices Af" are linearly independent. This ensures 
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the existence of a matrix close to the identity matrix such that a = tM^ for some 
unknown scalar quantity r. We may then apply Lemma 13.11 This concludes the proof 
of the theorem. □ 

We presented the above result for 7 in the vicinity of Idn in order to obtain a simple 
proof of a construction that satisfies the hypotheses of Lemma 13.11 and because the 
construction also appears in a later section. In fact, the result may be generalized as 
follows. 



Theorem 4.2 The results of Theorem 4-i hold for X the space of c sufficiently close to 



(ao,0,0), where ao is an arbitrary constant symmetric matrix satisfying (E]). 

Proof. The proof is very similar to that of the preceding theorem and is in some 
sense included in the proof of Theorem 14. 71 below, to which we refer for the details. The 
construction of Mi = 1 and Vj = Xj is the same as that of Theorem 14.11 The solutions 
Vn+m are then constructed as 

Vn+m = ^QmX 'X, 1 < TU < ^n{n + 1) - 1, 

with forming a family of M„ = ^n{n+ 1) — 1 linearly independent matrices that are 
orthogonal to Oq, or in other words, such that ao : Qij = 0. The linear combinations 6*™" 
are then constructed as in Theorem 14.11 with the matrices M"^ = Qm since V^'^vj = 
for 1 < j < n. This allows us to verify the hypotheses of Lemma [XT] globally on X = Xq 
for boundary conditions equal to the traces of the polynomials 1, Xj, \QmX ■ x, and by 
continuity for an open set of boundary conditions and for all coefficients c sufficiently 
close to (ao, 0, 0). □ 



4.2 Global reconstructions close to isotropic tensor 

Let us generalize the above result by assuming that a is in the vicinity of -y{x)Idn where 
7 is a scalar real- valued (hence positive) diffusion coefficient. We still assume that b 
is in the vicinity of 0. Also, c is an arbitrary complex-valued potential so that ([1]) is 
uniquely solvable. Then we have the following result. 

Theorem 4.3 Let j{x) e i^t+4+^(X) and c(x) G H'^+'^+'{X) for e > with &, in (Q 
bounded. Let X be the space of c sufficiently close to {'y{x)ldn,0,c{x)) in the sense of 
Lemma [Ol with m = 0. Let In = ^n{n + 3). There there exists an open set of {fi)i<i<j^ 
(in any topology of sufficiently smooth functions on dX ) such that 9Jtf from X to 2)^" is 
injective. Moreover, the stability result (120|) holds for Xq = X. 

Proof. The proof is based on the construction of complex geometrical optics solutions 
of the form 

u{x;p) = ^=e^-^{l + ij,{x)), (22) 

with p a complex- valued vector such that p ■ p = 0. We know that for 7 and c with 
the aforementioned regularity and for \p\ sufficiently large, then i/jp is of order \p\~^ in 

c\x) m- 
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In the construction above Lemma 13. we need to consider derivatives of ratios of 
solutions. We find that 

4^v4^=p-p+^, 4^v^^4^ = (p-pr+0, (23) 

u{X] p) u[x] p) u{X] p) u[x; p) 

with the vector tp bounded independent of (p, p) and matrix of order max(|p|, |p|) 
uniformly in x G X. 

Let us define pij = i{ei + iej) and define Uij = u{-; pij) as well as Uij = u{-; p*^). Note 
that Uij is asymptotically close to u*j as |€| — )■ oo but since c{x) may be complex valued, 
is not necessarily equal to u*j. 

We also define pij = epij as well as pi = e^pi2 and Ui = u{-; pi), with e'^t sufficiently 
large that contributions such as Lp and above remain negligible for the forthcoming 
constructions but e sufficiently small that p = pi or p = pij in fl23l) is so small that it 
does not modify the independence of the matrices M™ constructed below. 

Let us define Vj as follows 

v, = ^^^, ,^ = <lPl^^ ,>2. 

Ui Ui 

Since all solutions vj do not vanish for t sufficiently large, it is clear that (Vvi , . . . , Vvn) 
form a basis with Hij = Vvi ■ Vvj a matrix with a uniformly bounded inverse for x & X 
(with a bound that depends on i and e). Moreover, we find that 

Here and below, we denote by ~ equalities up to terms such as Lp and above that are 
asymptotically negligible as |t| — )■ oo as well as terms that are lower order in e. 
Now for 1 < 2 < j < n, we define 

Vij = — , SO that — Vvij ~ Pij, — Vij ~ pfj . 

Ui Vij Vij 

For l<z<n — 1, we construct 

V, = so that ^Vv, ~ p*,+„ ^V^'v, ~ iPh-^if- 

Each of the vectors "^Vij and V5j can uniquely be written in terms of the vectors Vvj. 
Let us define pi = p^2 pj = Pj-ij- Note that {Pj)i<j<n form a basis of C". Let us 
then introduce 

-Pij = e&^pk, -Pj,j+i = ^^ih- 

Here, the summation is over the index k. We find that all coefficients are of order 
e~^. Then we find that 

-Vv,, = &^Vv,, -Vv, = eiVv, for G^' ~ ^^QI 6^ ~ 

Now for these choices, we find that 

M^^ := + Q^^V^\, = v,,[pf - Q^^e^pf) ~ ^,,pg^ 
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Similarly, we have 
Note that 

jip- = i{ei^ej+ej^ei) + {ei®ei-ej^ej), ^'^^^ = -z(ej(8)ej+ej®ej) + (ej®ej-ej®ej). 

Therefore the matrices M^^ and constructed above are indeed linearly independent 
and as in the proof of Theorem 14.11 span a subspace of the vector space of symmetric 
matrices of dimension ^n{n — l)+n — l = |n(n + l) — 1. The above result obtained in the 
limit £ — > still holds for e sufficiently small. Moreover, once e is fixed, several terms 
of the form (p and f above become negligible when \t\ is sufficiently large. Therefore, 
for \t\ sufficiently large, there exists an open set of boundary conditions (fi) such that 
all the hypotheses of Lemma [3.11 are satisfied for all x & X. 

This proves the result for c of the form {'-f{x)ldn,0,c{x)). Now by continuity and 
Lemma 13.31 the same boundary conditions can be used to satisfy the requirements of 
Lemma 13.11 for all c sufficiently close to c. □ 



4.3 Global reconstructions in two dimensions 

The above reconstruction procedure has been proved to hold in the vicinity of 0, 0) 
or {'y{x)ldn,0,c{x)). In this section, we generalize the result to proving that global 
reconstructions are possible for coefficients in the vicinity of (7(x), 0, c(x)) where 7 is 
an arbitrary real-valued second-order elliptic tensor in dimension n = 2 and c{x) is a 
complex-valued potential: 

Theorem 4.4 Let X be the space of t sufficiently close to (7, 0, c) in the sense of Lemma 
rOI with m = with 7 of class H^^%X) and c of class H^^%X) for e> 0. Let I2 = 5. 
Then there exists an open set of boundary conditions f = {fi)i such that 9Jlf from X to 
2)^^ is injective. Moreover, the stability result (120|) holds for Xq = X. 

Proof. To simplify the notation, we set c = and leave the details to the reader to 
consider the case c 7^ as was done in the proof of Theorem 14. 3[ 

We prove that we can apply Lemma [XT] for Xq = X and (a, 6, c) = (7,0,0) for an 
open set of boundary conditions f. Then by continuity, the hypotheses of Lemma 13.11 
still hold for c sufficiently close to (7,0,0). As we did in the proof of the preceding 
theorem, we can assume that 7 is smooth since by an application of Lemma 13. 3[ the 
result can be extended to any 7 satisfying the regularity hypotheses of Lemma [3. 31 with 
m = 0. 

This global reconstruction works only in two dimensions of space and for real valued 
tensors 7. The reason is that global complex geometrical optics solutions can be con- 
structed in two dimensions by means of appropriate quasiconformal maps. Such results 
do not hold in general in dimension n > 3. 

Let Go be the identity conformal structure and G = G{z) the conformal structure 
given by 7. Then there is a diffeomorphism (p from C — !■ C, unique after normalization 
at infinity, such that [5] 

D(t)\z)J-\z,(f){z))D(j){z) = G{z). 
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Then with (f)*u = u o (p we have that 

V ■ J-^V(j)*u = V ■ GVu o (j). 

Identifying z = xi + 1x2 and ) and taking (t)*u{x) = \/j{(f){x))e'^'^ , we find 

that u(x) = \/^{x)e'^''^^^^ with ^{x) = {ipi(x), (p2{x)) a diffeomorphism with ip~^ = 0. 
In other words, we construct CGO solutions for anisotropic media of the form 

Up{x) = yj(x)e''-'^(")(l + M^))- 

Let pi = k{iei + 62) and p2 = k{ie2 — Ci). Let Uj = Up., j = 1,2. As in the preceding 
section, we also define Uj = Up* for j = 1,2. Then we find 

k^^Vui = {iVipi + 'V(p2)ui + Ci, k~^Vu2 = i{iV(pi + V(p2)u2 + (2- 
We thus find 

k~\2'^Vu2 = ik^^u^^Vui + 0(1) = -Vv?! + zVv22 + 0(1)- 

Here, we are decomposing Vm2 over Vui and V-Ui, which form a basis for k sufficiently 
large. Then 

M = ik~'^u{^V ® Vui - k-'^U2^V ® Vu2 

~ (i + 1) (V^f - V^T) + (V(/7i ® Vip2 + V(/?2 ® V(/?i) , 

in the limit /c — t- 00. In the same way that we have decomposed Vn2 over Vui and 
Vui above, we can decompose Vu2 over Vtti and Vui as well. In the limit A; — 00, the 
matrix M* will thus be given by the complex conjugation of the above matrix. This 
proves that by in the limit — )■ 00, the matrices M"* that we construct are given by the 
real and imaginary parts of M: 

M± = (Vv^f - V^f) ± (V(/Ji ® V^2 + Vip2 ® Vv?i) . 

After change of coordinates, we obtain the two matrices: 

Ml = V<^f ^ - Vyjf ^ M2 = V<^i ® V<^2 + V(/72 ® Vy?! 

which we want to be non trivial and linearly independent. The above matrices Mi 2 are 
those obtained in the limit k — )■ 00. This means that for k sufficiently large, the two 
constructed matrices Mi^2 from ui and U2 will be close to their limits and hence satisfy 
the same properties of linear independence. 

Now we observe that a = V(fi and b = V(p2 are linearly independent since (f is a. 
diffeomorphism. And a® a, b^b, a®b + b®a are basis elements for symmetric matrices. 
Thus Ml has coordinates (—1, 1,0) while M2 has coordinates (0,0, 1) in that basis. As 
a consequence, both matrices Mi and M2 are linearly independent, in the limit A; — > 00 
as well as for k sufficiently large. Note that the independence is uniform in x G X for 
k sufficiently large. 

This shows that the hypotheses of Lemma 13.11 are satisfied for Xq = X. Such a 
calculation holds for any set of coefficients close to (7, 0, 0). A very similar proof applies 
to c in the vicinity of (7, 0, c) as stated in the theorem. This proves the theorem. □ 
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4.4 Global reconstructions with redundant measurements 

In this section, we show that reconstructions are possible for essentially arbitrary (suf- 
ficiently smooth) coefficients c. However, the construction of the matrices M™' becomes 
local. We thus need to use a number of internal functionals / that is potentially much 
larger than /„, although we do not expect this large number of coefficients to be neces- 
sary in practical inversions. 

The local constructions require that certain properties of linear independence be sat- 
isfied. Such conditions will be satisfied for well-chosen illuminations fj on the boundary 
dX. The control of the linear independence from the boundary is obtained by means 
of a Runge approximation; see Lemma 14.81 below. This step requires that the operator 
L = V • aV + b-'V + c satisfy a unique continuation principle, which we state as follows: 

Property 4.5 (Unique Continuation) We say that L satisfies the unique continu- 
ation principle when Lu = on X\Xq with u = on dX and n ■ Vm = on dX 
implies that u = on X\Xo, where Xq is an arbitrary sufficiently smooth open domain 
Xo CC X. 

For unique continuation results, we refer the reader to [H [16] and the theoretical results 
we shall use here [lOl Theorem 17.2.1]. The latter result states that L satisfies the unique 
continuation principle 14.51 when the principal symbol of L given by p(x,^) = a{x)^ ■ ^ is 
such that: 

(i) a{x) is Lipschitz continuous, 

(ii) For E M"\{0}, the quadratic equation p{x,^ + tN) = in the 
variable r G C admits a double root r if and only if ^ -|- tN = 0. 

Then we have the following lemma: 

Lemma 4.6 Let p{x, ^) = a{x)C, ■ C, be the principal symbol of L, which we assume is 
elliptic. 

In dimension n > 3, the quadratic equation p{x, ^ + rN) = for ^, N E ]R"\{0} never 
admits a double root r unless ^ + tN = 0. 

In dimension n = 2, the same result holds when in addition ([2]) is satisfied. 
In all these cases, L thus satisfies Property \4 .5\ when a is Lipschitz continuous. 

Proof. The proof is essentially given in [TUl Lemma 17.2.5]. In dimension n > 3, 
the equation p{x, ^ + tN) has one root with 53r > and one root with 5Rr > 0. In 
dimension n = 2, the equation p{x, C, + tN) has a double root at a fixed point x G X if 
and only if we have p{x,^) = {l{x) ■ ^)^ for some complex-valued vector I = Ir + ik E C" 
(with Ij. = JJZ and li = Ql) to preserve ellipticity (note that such quadratic forms cannot 
be elliptic in dimension n > 3). But then p{x,^) = {Ir ■ 0^ ~ i.h ■ + i^r ■ Ch' C so that 
the real part of a is not elliptic. This proves the lemma. □ 

With this result, we can now state the main theorem of the paper. 

Theorem 4.7 Let X be the space of coefficients c such that (6, c, V ■ a) are of class 
C'''"(X), a is of class C^'^{X), and such that ([2]) holds. Then there exists I > In and 
an open set (for the topology of C'^'°'{dX)) of boundary conditions f = (/j)i<i</ such 
that OJtf from X to 2)'^ is injective. Moreover, the stability result ( l20i) holds for Xq = X . 



14 



Proof. We decompose the proof into three steps: we first construct local solutions 
assuming that the coefficients are constant. We then extend the local constructions 
to the case of non-constant coefficients. We finally apply the Runge approximation to 
obtain an open set of boundary conditions such that the hypotheses of Lemma 13.11 are 
satisfied locally. Local constructions are then patched together to provide global stable 
and unique reconstructions. 



Problem with constant coefficients. Let first Xq be a point inside X, which by 
change of coordinates we call 0. Let us define {ao,bo,Co) = (a(0),6(0) + V ■ a(0),c(0)). 
We then look for solutions of the constant coefficient equation 

LqU := ao : V®^u + feo ■ Vu + cqU = 0. (24) 

We look for solutions approximately of the form 

p = ^Qx ■x + p-x + d, Vp = Qx + p, V®^p = Q. 

In order for p to satisfy the equation at x = 0, we need to find {Q, p, d) such that 

ao : Q + bo ■ p + Cod = 

We construct /„ = 1 + tt, + ln(n + 1) — 1 such solutions below. We then realize that 
LqP = 0{x). Let To be sufficiently small and let us define 

LoU = in 5(0, ro), u = p onas(0,ro). (25) 

For To sufficiently small, the derivatives up to order two of u and p are very close. 
The linear independence of the structures constructed below with the polynomials p at 
x = therefore still holds for the corresponding structures constructed with the elliptic 
solutions u in ( l25l) . 

We call the first solution Uq obtained by defining 

rf=l, p = 0, Q- 



Note that Uq does not vanish in a sufficiently small neighborhood of (and can be 
normalized so that Uo(0) = 1). We next define the solutions Uj for 1 < j < n. The 
vector ho = bor + iboi is after a rotation if necessary in the span of ei and 62- We thus 
write bo = pei + 1^62 for p and p in C. For j = 1, 2, we define 



For J > 3, we define 

d = 0, p = Cj, Q = 0. 

The solutions Uj are therefore constructed such that Vpj = ej at x = 0. Moreover, we 
find that Vr^ = at a; = as well since Vpn = at x = 0. We thus obtain that V— 
form a basis of in a sufficiently small neighborhood of 0. 
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Finally, for n + 1 < j < n+ + 1) — 1, we define d = and p = and choose the 
matrices Qj such that they form a free family of symmetric matrices that are orthogonal 
to ttg, the complex conjugate of uq. This free family has dimension ^n{n + 1) — 1. This 
implies that (oq : Qj) = 0. 

In the construction above Lemma 13.11 it is V^^f j+m for vj = Uq^Uj that is used to 
form a free family of dimension |n(n + 1) — 1. We verify that 

V®2^ = -V^\+m + Vu,+„ ® V- + V- ® Vu,-_,„ + u,+,„V^2-. 

Uo Uo Uo Uo Uo 

We verify that both sides equal Qj at x = when Uj is replaced by pj. 

Let now OJ^ and M™ be defined as above Lemma 13.11 We verify that M"^ is close to 
Qj+m at X = 0. By continuity, the matrices M™' are therefore linearly independent in a 
ball Xq = B{0, To) for tq > sufficiently small. This shows that on that ball, the family 
of matrices as constructed above Lemma 13.11 satisfy the hypotheses of that lemma. 
All other hypotheses of that Lemma are therefore satisfied for the family Uj. 

Problem with non-constant coefficients. We now return to the full problem and 
look for solutions of the form: 

a{x):V^'^u+{b + V ■a)-Vu + cu = in Xq, u = u on dXo, (26) 

where Xq = B{0, ro) is a ball whose radius ro is equal to or smaller than the value chosen 
in the construction of u. Let w = u — u. We find 

ao : V^^w + bo ■ Vw + CqW = (a - Oq) : V®^ii + (6 + V ■ a - 6o) ■ Vn + (c - Cq)u 

on Xq with w; = on dX^. By assumption on the coefficients and -u, the above right- 
hand side is bounded uniformly by Tq on the ball Xq. We deduce from elliptic regularity 
results for complex-valued coefficients [HI Chapter 6] that 

lkllc2>-(Xo) < Cr'^, 

for some positive constant C. Thus for ro sufficiently small, we find that the functions 
Uj defined as solutions for fl2^ with boundary conditions Uj are arbitrarily close to Uj 
in the sense for ro sufficiently small. This proves that Uq remains non-vanishing and 
close to 1 on Xo, that for 1 < j < n remain linearly independent, and that the 
matrices M™ constructed above Lemma 13.11 satisfy the independence properties stated 
in that lemma. 

Continuation to the boundary. So far, we have constructed solutions u that are 
defined on Xq = B{0, ro). We need to construct solutions on the whole domain X such 
that their restrictions on X^ = B{0, /iro) is a sufficiently accurate approximation of u 
for < fi < 1. We need the following Runge approximation property, following |T5]; see 
also m\. 

Lemma 4.8 (Runge approximation) Let L be an operator satisfying the unique con- 
tinuation property of Cauchy data on X as described above. 
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Let uq be a solution of Luq = on Xq and let = B{0, ^tq) for < /i < 1. Then 
for each e > 0, there is a function fe E H^(dX) such that the solution of Lue = on 
X with Us = fe on dX is such that 

\\ue- uo\\c2,c(^x^) < e. (27) 

Proof [Runge Lemma]. Let E = {u E H^{Xq), Lu = in Xq} and F = {u^Xo^ ^ ^ 
H^{X), Lu = in X} be linear subspaces of L'^{Xq). We wish to prove that F = E 
for the strong topology. By Hahn Banach, this means that for all / G L^{Xo), then 
(/, u) = for all M G F implies that (/, u) = for all u E E. 

Let us extend / by outside Xq and still call / the extension on X. Define then 

L*v = f in X, v = on dX. 

Here L* = V ■ a*V — 6* ■ V + c* is the formal adjoint to L. Note that v is well-defined 
since ©c in ([6]) is assumed to be bounded. Integrations by parts show that 

{Lu, v) — {u, L*v) = / (an ■ Wuv* — an ■ Vv*u — b ■ nuv*)da = an ■ Wv*uda. 

JdX JdX 

Since this holds for any function u E H^{X) and hence for any u\dx E H^(dX), we 
deduce that a*n ■ Vv = on dX. We thus find 

L*v = in X\Xq, V = and a*n ■ Vv = on OXq. 

We use the unique continuation assumption to deduce that v = in {X\Xq) so that 
V = in the H 2(^8X0) sense and a*n ■ Vv = in the H~2{dXQ) sense. For any u E E, 
we thus find that (/, u) = 0, which thus proves that F = E. This shows that u — Uq is 
arbitrarily small in L^(Xo). Now regularity results as they are written for instance in 
p31 Theorem 6.2.5] for elliptic problems with complex coefficients such that (a, b, c, V-a) 
are of class C°'" for a > (see also [TOl Theorem 17.2.7]), allow us to conclude that (1271) 
holds. Indeed, we have an equation L{uo — u^) = on Xq. We then get the required 
interior regularity of Uq — u^ in C^'"(i?(0, ^rg)) for all ^ <1. □ 

We now conclude the proof of Theorem 14.71 The uniqueness to the Cauchy problem 
is guaranteed by Lemma 14.61 We have obtained, using the Runge approximation, the 
construction of a family Uj for an open set of boundary conditions fj such that the 
hypotheses of Lemma 13. II are satisfied on Xq. It remains to cover X by a finite number 
of balls of radius /xro (for /i < 1 as necessary to apply the Runge approximation result) 
and to apply Lemma 13.11 globally on X and obtain a unique and stable reconstruction 
of z on X. □ 



5 Applications to coupled-physics inverse problems 

The salient feature of coupled-physics inverse problems (also known as hybrid inverse 
problems) is that they involve a high resolution modality and a high contrast modality 
to obtain a coupled (hybrid) modality imaging combining both high contrast with high 
resolution. We consider three such families of coupled physics inverse problems that 
may be modeled by the theory developed in the preceding sections; quantitative photo- 
acoustic tomography (QPAT), transient elastography (TE), and the mathematically 
similar modality called magnetic resonance elastography (MRE); see |1] for a review on 
hybrid inverse problems. 
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5.1 Quantitative Photo- Acoustic Tomography 

The first modality we consider is called quantitative photo-acoustic tomography (QPAT). 
The high contrast modality is optical tomography. The ultimate objective of QPAT is 
the reconstruction of the optical coefficients in an elliptic equation. 
Radiation propagation is modeled by the following equation 

- V ■ 'jVuj + auj = in X, uj = fj on dX. (28) 

Here, 7 is the real- valued diffusion tensor and a the real- valued absorption coefficient. 

The high resolution modality is ultrasound. A first well posed inverse wave (ultra- 
sound) problem is solved to reconstruct internal functionals of the unknown coefficients. 
This first step of QPAT provides access to the following internal functionals O [6l [71 [19] 

Hj{x) = r{x)a{x)uj{x) inX (29) 

Here T{x) is the Griineisen coefficient, which is assumed to be known in this paper and, 
therefore, without loss of generality assumed to equal 1. We assume that 1 < j < /, 
with / the number considered in the preceding sections. We assume that all coefficients 
are known on dX and that /i > on dX so that i^i > by the maximum principle. 
Then multiplying the above equation for ui by uj and for uj by Ui and subtracting the 
results, we get 

-V ■ hui)V^ = -V ■ hul)V§- = 0. 
ui Hi 

Therefore if the / — 1 conditions (/2, • • • , //) boundary conditions are chosen as in the 
preceding section, we obtain that ju^ can be uniquely and stably reconstructed. Indeed, 
we are here in the setting where 6 = 0, which allows one to reconstruct the gauge and 
hence the whole diffusion tensor 7 as indicated in Corollary 12.21 
Now the equation for ui may be recast as 

- V ■ (7M?)V— = Hi in X, — = ^ on dX. (30) 

Ui Ui fi 

This uniquely determines Ui and hence 71 in a stable fashion. Since Hi = aui, this 
also determines a uniquely and stably. This concludes the derivation of the unique and 
stable reconstruction of (7, a) from QPAT measurements when the Griineisen coefficient 
is known. 

Note that the same elliptic equation ( l28l) with cr = has been used to reconstruct a 
scalar diffusion coefficient from knowledge of u by solving the transport equation ( l28l) 
for 7 with applications in underground water flows [H [IT] . 

5.2 Coupled-physics methods based on Elastography 

In this section, the high contrast modality is elastography; see [12] and reference there 
for more details. The elastic (stiffness) properties of tissues are to be reconstructed. 
We assume here that the elastic displacements are modeled by a time-harmonic scalar 
equation of the form 

V ■ 'j{x)\/uj + uj^p{x)uj = in X, uj = fj on dX. (31) 
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Here, 7 is a tensor- valued, possibly complex- valued. Lame parameter and p is a density 
that may also be complex-valued in full generality. The reason for these coefficients 
to be complex-valued is that elastic waves are attenuated by various dispersion effects. 
In the frequency domain, such attenuation effects take the form of complex-valued co- 
efficients. Elastographic tomography was one of the main motivations to consider the 
reconstruction of complex-valued coefficients in the preceding section. 

In transient elastography (TE), the high resolution modality is again ultrasound. 
As comparatively slow elastic waves propagate through the domain X of interest, ul- 
trasound measurements are used to infer the internal displacements, i.e., the solution 
of the Helmholtz equation ( 1311) . In Magnetic resonance elastography (MRE), the high 
resolution modality is magnetic resonance. Along with elastic displacements, proton 
displacements occur that can be measured by an MRI machinery. 

In both modalities, the internal functionals obtained by ultrasound in TE and by 
magnetic resonance imaging in MRE are given by the displacements: 



This is exactly the setting considered in Theorem 12.11 Note, however, that in many 
applications of elastography, the scalar model considered here is not sufficiently accurate. 
Generalizations to more precise models of linear or nonlinear elasticity then need to be 
developed. 
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